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C . Abstract. We prove global existence of nonnegative weak solutions to a degenerate parabolic 

1 ^ ' system which models the interaction of two thin fluid films in a porous medium. Furthermore, we 



show that these weak solutions converge at an exponential rate towards flat equilibria. 

1. Introduction 

In this paper we consider the following system of degenerate parabolic equations 



o 

(N 

< 

j dtf= d.Afd.f) + RdAfd.h), ^ ,n ^ r^ 

{ G (0,00 X (0, L), (1.1) 

\dth= d,ifdJ) + R^d,[{h-f)d,h] + RdAfd.h), 

which models two-phase flows in porous media under the assumption that the thickness of the two 
^ ' fluid layers is small. Indeed, the system (1.1) has been obtained in [4] by passing to the limit of 

. small layer thickness in the Muskat problem studied in [3] (with homogeneous Neumann boundary 

condition). Similar methods to those presented in [4] have been used in [6] and [8], where it is 
■ rigorously shown that, in the absence of gravity, appropriate scaled classical solutions of the Stokes' 

and one-phase Hele-Shaw problems with surface tension converge to solutions of thin film equations 

dth + d^{h^d^h) = 0, 

with a = 3 for Stoke's problem and a = 1 for the Hele-Shaw problem. 

In our setting / is a nonnegative function expressing the height of the interface between the fluids 
^ I while h> f is the height of the interface separating the fluid located in the upper part of the porous 

^ ' medium from air, cf. Figure 1. We assume that the bottom of the porous medium, which is located 

at y = 0, is impermeable and that the air is at constant pressure normalised to be zero. The 
parameters R and R^ are given by 

R := — — — and i?^ := ^^R , 

P- - P+ P+ 

where /i_ [resp. ^j^] denote the density and viscosity of the fluid located below [resp. above] 
in the porous medium. Of course, we have to supplement system (1.1) with initial conditions 

/(0) = /o, MO) = /io, xG(0,L), (1.2) 
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Figure 1. The physical setting 

and we impose no- flux boundary conditions at x = and x = L: 

dj = d^h = 0, x = 0,L. (1.3) 

It turns out that the system is parabolic if we assume that /ig > /o > and R > 0, that is, the denser 
fluid lies beneath. Existence and uniqueness of classical solutions to (1.1) have been established in 
this parabolic setting in [4]. Furthermore, it is also shown that the steady states of (1.1) are flat 
and that they attract at an exponential rate in solutions which are initially close by. 

In this paper we are interested in the degenerate case which appears when we allow /o = and 
ho = fo on some subset of (0, L). Owing to the loss of uniform parabolicity, existence of classical 
solutions can no longer be established by using parabolic theory and we have to work within an 
appropriate weak setting. Furthermore, the system is quasilinear and, as a further difficulty, each 
equation contains highest order derivatives of both unknowns / and h, i.e. it is strongly coupled. 
In order to study the problem (1.1) we shall employ some of the methods used in [2] to investigate 
the spreading of insoluble surfactant. However, in our case the situation is more involved since we 
have two sources of degeneracy, namely when / and g := h — f become zero. It turns out that by 
choosing {f,g) as unknowns, the system (1.1) is more symmetric: 

j dtf= {l + R)dAfd,f) + RdAfd.g), ,,,,,, 

< {t,x) G (0,oo) X (0,L), (1.4) 

[ dtg = R^dx (gdxf) + R^^d^ (gd^g) , 

since, up to multiplicative constants, the first equation can be obtained from the second by simply 
interchanging / and g. Corresponding to (1.4) we introduce the following energy functionals: 



Siif,9) ■■= t 
Jo 



(/ln/-/ + l) + ^(<7ln5-5 + l) 



dx 



and 

dx. 



£2{f,g):= [ [f + R{f + g? 
Jo 



It is not difficult to see that both energy functionals £i and £2 dissipate along classical solutions of 
(1.4). While in the classical setting the functional £2 plays an important role in the study of the 
stability properties of equilibria [4], in the weak setting we strongly rely on the weaker energy £1 
which, nevertheless, provides us with suitable estimates for solutions of a regularised problem and 
enables us to pass to the limit to obtain weak solutions. Note also that £1 appears quite natural 
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in the context of (1.4), while, when considering (1.1), one would not expect to have an energy 
functional of this form. 

Our main results read as follows: 

Theorem 1.1. Assume that R > 0, > 0. Given fo,go G L2((0, L)) with fo>0 and go > there 
exists a global weak solution {f,g) of (1.4) satisfying 

(i) f>0,g>Otn (0,r) x (0,L), 

(ii) f, 9 G ^oo((0, T), L2((0, L))) n L2((0, T), H\{0, L))), 
for all T > and 



fmdx 



g{T)^ dx 



foipdx 



goipdx 



T 

Jo 



-R. 



((1 + R)fdJ + Rfd^g) dx dt, 

T f-L 

(gdxf + gdxg) d^ip dx dt 



JO 



for all ip S W^{{0, L)). Moreover, the weak solutions satisfy 
(a) ||/(T)||i = ||/o||i,||5(T)||i = |bo||i, 



(6) £i{f{T),g{T)) + 

(c) f2(/(T),5(T)) + 
for almost all T G (0,oo). 



T 

JO 

T 



dxdt < £i{fo,go), 



JO 



f ((1 + R)dJ + Rd^gf + RR^g{dJ + d^gf dxdt < £2{fo,9o) 



Remark 1.2. If /q = for instance, a solution to (1.4) is {0,g) where g solves the classical porous 
medium equation dtg = R^dx {gdxg) in (0, oo) x (0, L) with homogeneous Neumann boundary 
conditions and initial condition g^. 

Additionally to the existence result, we show that the weak solutions constructed in Theorem 1.1 
converge at an exponential rate towards the unique flat equilibrium (which is determined by mass 
conservation) in the L2— norm: 

Theorem 1.3 (Exponential stability). Under the assumptions of Theorem 1.1, there exist positive 
constants M and uj such that 

2 , „r, 2 



fit) 



L 



fodx 







+ 



git) 



godx 







< Me 



for a.e. t > 0. 



Remark 1.4. Theorem 1.3 suggests that degenerate solutions become classical after evolving over 
a certain finite period of time, and therefore would converge in the i/^— norm towards the corre- 
sponding equilibrium, cf. [4]. 

The outline of the paper is as follows. In Section 2 we regularise the system (1.4) and prove that 
the regularised system has global classical solutions, the global existence being a consequence of 
their boundedness away from zero and in H^{{0, L)). The purpose of the regularisation is twofold: 
on the one hand, the regularised system is expected to be uniformly parabolic and this is achieved 
by modifying (1.4) and the initial data such that the comparison principle applied to each equation 
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separately guarantees that / > e and g > e for some e > 0. On the other hand, the regularised 
system is expected to be weakly coupled, a property which is satisfied by a suitable mollification of 
dxQ in the first equation of (1-4) and d^f in the second equation of (1.4). The energy functional £i 
turns out to provide useful estimates for the regularised system as well. In Section 3 we show that 
the classical global solutions of the regularised problem converge, in appropriate norms, towards a 
weak solution of (1.4), and that they satisfy similar energy inequalities as the classical solutions of 
(1.4) for both energy functionals £i and £2- Finally, we give a detailed proof of Theorem 1.3. 

Throughout the paper, we set Lp := Lp{{0,L)) and := W^{{0,L)) for p G [l,oo], and 
:= //^"((O, L)) for a G [0, 1]. We also denote positive constants that may vary from line to line 
and depend only on L, R, R^, and {fo,go) by C or Q, i > 1. The dependence of such constants 
upon additional parameters will be indicated explicitly. 

2. The regularised system 

In this section we introduce a regularised system which possesses global solutions provided they 
are bounded in and also bounded away from zero. In Section 3 we show that these solutions 
converge towards weak solutions of (1.4). 

We fix two nonnegative functions /o and go in L2 (the initial data of system (1.4)) and first 
introduce the space 

Hi := {/ G F2((0,L)) : dJiO) = dJiL) = 0}. 

We note that for each e > 0, the elliptic operator (1 — e^9^) : — )• L2 is an isomorphism. This 
property is preserved when considering the restriction 

(1 - e^dl) : {/ G C2+-([0,L]) : dj{0) = dJ{L) = 0} ^ C^{[0,L]), a G (0, 1). 

Given f,g&L2,we let then 

:= (1 - e'dlr'f, := (1 - e'dlr'g, (2.1) 

and consider the following regularised problem 

f dtfe= {l + R)dAfedJe)+RdAife-e)d,Ge), ,,,,,, , , , , 

{ (t,x) G (0,00) X (0,L), (2.2a) 

( dtge= Rf^dx {{ge - £)dxFe) + R^,d.x {gedxge) , 

supplemented with homogeneous Neumann boundary conditions 

dxfe = dxge = x = 0,L, (2.2b) 

and with regularised initial data 

fe{0) = he ■■= (1 - e'52)-Vo + £, 9e{0) = 50e := (1 - e^d^y^go + e. (2.2c) 

Note that the regularised initial data (/oe,5'Oe) £ x and, invoking the elliptic maximum 
principle, we have 

foe > e, goe > £■ (2.3) 
Letting Fq^ := (1 — e^5^)"^/o and Goe := (1 — £'^dx)~^go, we obtain by multiplying the relation 
i^Oe — e^(9^Foe = /o by and integrating over (0, L) the following relation 



[ foFoedx<\\fo\\2\\Foe\\2, 

Jo 
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which gives a uniform L2-bound for the regularised initial data: 

Wfoeh < II/0II2 + eVZ and Wgoeh < Ibolb + eVZ. (2.4) 

Concerning the solvability of problem (2.2), we use quasilinear parabolic theory, as presented in 
[1], to prove the following result: 

Theorem 2.1. For each e € (0, 1) problem (2.2) possesses a unique global nonnegative solution 
Xe ■■= [fe^ge) with 

fe,ge G C([0,cx)),i7i)nC((0,oo),Fi)nCi((0,oo),L2). 

Moreover, we have 

fe>£, > £ for all (i, x) G (0, 00) X (0, L). 

In order to prove this global result, we establish the following lemma which gives a criterion for 
global existence of classical solutions of (2.2): 

Lemma 2.2. Given e G (0, 1), the problem (2.2) possesses a unique maximal strong solution = 
ife,ge) on a maximal interval [0,r+(e)) satisfying 

fe,ge€C{[0,T+{e)),H^)nCmT+{e)),Hl)nC\{0,T+{e)),L2). 
Moreover, if for every T < T^{£) there exists C{£,T) > such that 

fe>e/2 + C{£,Tr\ ge>e/2 + C{e,T)-\ and max ||X,(t)||^i < C(e,T), (2.5) 

te[o,T] 

then the solution is globally defined, i.e. Tj^[e) = 00. 

Proof. Let £ G (0, 1) be fixed. To lighten our notation we omit the subscript e in the remainder 
of this proof. Note first that problem (2.2) has a quasilinear structure, in the sense that (2.2) is 
equivalent to the system of equations: 

{ dtX + A{X)X = F{X) in (0,00) X (0,L), 

BX = Q on (0,00) X {0,L}, (2.6) 

X(0) = Xo on (0,L), 

where the new variable is X := (/, g) with Xq = (/o, 50)1 and the operators B and F are respectively 
given by 

■ R{f - e)d.,G \ 



BX := d,X, F{X) := 
Letting 

a{X) := 



R^{g-e)d^F 



{l + R)f \ 
R^g )' 

the operator A is defined by the relation A{X)Y := —dx{a{X)Y). We shall prove first that (2.6) 
has a weak solution defined on a maximal time interval, for which we have a weak criterion for 
global existence. We then improve in successive steps the regularity of the solution to show that 
it is actually a classical solution, so that this criterion guarantees also global existence of classical 
solutions. 
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Given a E [0, 1], the complex interpolation space := [L2>-f^g]a is known to satisfy, cf. [1], 

^2. ^ f ' " < 3/4, 

^ \ {/ G : d^f = for X = 0, L} , a> 3/4. 

Furthermore, for each /? G (1/2, 2] we define the set 

■■= {f^H§: e/2 < /}, 

which is open in H^. Choose now 7 := 1/2 — 2^ > 0, where ^ G (0, 1/18). We infer from (2.2c), that 
^0 £ ^ • III order to obtain existence of a unique weak solution of (2.6) we verify the 
assumptions of [1, Theorem 13.1]. With the notation of [1, Theorem 13.1] we define 

{a, s, r, r) := (3/2 - 3^, 1 + ^, 1 - -0, 2a := 3/2 - 3^ 

Since 1 — ^ — 1/2 > 7, we conclude that Vg~^ C C^{[0,L]), meaning that the elements of a{X) 
belong to C^{[0,L]) for X G Vg"^ x V^^^. Since a{X) is positive definite, we conclude in virtue of 

7 > 2S - 1 that 

(Afi)GCi"(V^-«xV^-«,£:"(0,L)), 

the notation £°'(0,L) being defined in [1, Sections 4 & 8]. Moreover, since {1 — e^d^)^^ G C{L2,H^) 
we also have 

F G C^^iV^^^ X V^~^,H~^), 

whereby, in the notation of [1], = H^^ because |^| < 1/2, cf. [1, eq. (7.5)]. Thus, we find all 
the assumptions of [1, Theorem 13.1] fulfilled, and conclude that, for each Xq G Vg, there exists a 

unique maximal weak Hq —solution X = {f,g) of (2.6), that is 

/,ffGC([o,r+),v^-«)nC((o,r+),i/^/'-«)nci((o,r+),ifg^/'-3«), 

and the first equation of (2.6) is satisfied for all t G (0,T+) when testing with functions belonging 
to h'J'+'^. m oreover, if X|[0,T] is bounded in x Hq and bounded away from aVg"^ for all 
T > 0, then T+ = 00, which yields the desired criterion (2.5). 

We show now that this weak solution has even more regularity, to conclude in the end that 
the existence time of the strong solution of (2.6) coincides with that of the weak solution (of 
course they are identical on each interval where they are defined). Indeed, given 5 > 0, it holds 
that X G C([5,r+),i/g/^"^^) n C^{[6,T+),H~^^^~^^). Hence, if < 2/)< 2, we conclude from [1, 
Theorem 7.2] and [7, Proposition 1.1.5] that X is actually Holder continuous 

X G CP{[6,T+),hI^^^^^^^''). 

Choosing p := ^ and := 1 — 6^ > 0, we have that H'^^'^ '^^ ^ ^"^([0, -^^j), so that the elements 
of the matrix a{X{t)) belong to C^([0, L]) for all t G [5, r+). Defining 2/2 := 2 - 8^ > 0, we observe 
that /i > 2/i — 1 and 

(^(x),fi)GC''([5,r+),f^((o,L))). 

Finally, with 1v := 3/2 + our choice for ^ implies 

iX{5),F{X)) G Hl''-^ X CPi[6,T+),H^^) c H^''-^ x C''i[5,T+), H^''-^ 
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and the assertions of [1, Theorem 11.3] are all fulfihed. Whence, the hnear problem 

dtY + AiX)Y = F{X) in (,5,r+) X (0,L), 

BY = on (5,r+) X {0,L}, (2.7) 

Y{d) = X{6) on (0,L), 
possesses a unique strong iJ^'^— solution Y, that is 

Y G c([5,r+),/?2-2 ^ ^2.^2) n c{{5,T+),Hf X i/i^) n ci((<5,r+),//2.-2 ^ ^2.-2^^ 

In view of [1, Remark 11.1] we conclude that both X and Y are weak H^^ —solutions of (2.7), 
whence we infer from [1, Theorem 11.2] that X = Y, and so 

/,<7GC([5,r+),V^'+«)nC((5,r+),i7^/'+«)nC7i((<5,r+),i/g^/'+«). 

Interpolating as we did previously and taking into account that 5 was arbitrarily chosen, we have 
G C^([5,r+),Ci+^([0,L])) if we set AO < ^. Hence, we find that 

{X{6),{A{X),F{X))) e L2 X C\[6,T+),n{Hl x hI,L2 x L2) x (L2 x L2)), 

where ^{{H^ x ff^, L2 x L2) denotes the set of linear operators in L2 x L2 with domain x //g 
which are negative infinitesimal generators of analytic semigroups on L2 x L2, and, in virtue of [1, 
Theorem 10.1], 

X G C{{5,T+),Hl X Hl)r^C\{5,T+),L2 x L2) 
for all 5 G (0,^4.). Hence, the strong solution of (2.6), which is obtained by applying [1, Theorem 
12.1] to that particular system, exists on [0, T-|_) and the proof is complete. □ 

The remainder of this section is devoted to prove that T+(e) = 00 for the strong solution {fe,ge) 
of (2.2) constructed in Lemma 2.2. In view of Lemma 2.2, it suffices to prove that {fs,ge) are a 
priori bounded in and away from zero. Concerning the latter, we note that we may apply the 
parabolic maximum principle to each equation of the regularised system (2.2a) separately. Indeed, 
owing to the boundary conditions (2.2b), the constant function {t,x) 1— )• e solves the first equation 
of (2.2a) and (2.2b) while we have foe > e by (2.2c). Consequently, > e in [0,T+(e)) x [0,L] and, 
using a similar argument for i^e, we conclude that 

fs>e, ge>e for all (t, x) G [0, r+(e)) x (0, L). (2.8) 

Next, owing to (2.2b) and the nonnegativity of and gg, it readily follows from (2.2a) that the 
Li— norm of and is conserved in time, that is, 

\\fe{t)\\i = WkeWi = ll/olli + eL, \\ge{t)\\i = Ml = lbo||i + sL (2.9) 

for all t G [0, T4_(e)). The next step is to improve the previous L^— bound to an bound as 
required by Lemma 2.2. To this end, we shall use the energy £1 for the regularised problem, see 
(2.13) below. As a preliminary step, we collect some properties of the functions {F^,Ge) defined in 
(2.1) in the next lemma. 

Lemma 2.3. For all t G (0,T+(e)) 

\mt)\\2<\\fem2, ||G,(t)||2< ||9.(t)||2, (2.10) 
\\dMt)h < \\dxfe{t)\\2, \\d^Gem2 < ||5.5e(t) lb , (2.11) 

e \\dlF,{t)h < \mem2, e \\dlGS)\\2 < \\d,ge{t)h. (2.12) 
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Proof. The proof of (2.10) is similar to that of (2.4). We next multiply the equation F^ — e'^d'^F^ = 
by —d'^Fe, integrate over (0, L) and use the Cauchy-Schwarz inequality to estimate the right-hand 
side and obtain (2.11) and (2.12). □ 

Lemma 2.4. Given T G (0, T+(e)), we have that 



£i{MT),g,{T)) + £ £ + -^—\d,g^\^^ dxdt < £i{fM,9em. 



(2.13) 



Proof. Using (2.11) and Holder's inequality, we get 



j^Slifs, Qe) = £ dt Ue log(/.)) + f-^dt {ge log{g,)) dx 



(1 + R)\dM^ + R^^d^fAGe + R^^d^gAF, + R\d,gs\^ 



fe 



96 



dx 



<-(l + R)\\d,fe\\l + RWd.fehWd.Geh + R\\d,geh\\d,Fe\\2 - R\\d,g. 



'e\\2 



— 7:\\^xfe\\2 



1 + 2R, 



9x/e||2 ~ 2i?||(9a;/e||2||9a:5e||2 + -R||c^z9e||2 



< 



1, 



dxfi 



R 



211— ^£ll2 i^2R' 
Integrating with respect to time, we obtain the desired assertion. 



:\\dx9. 



2 

'e|l2- 



□ 



Since z In z — z + 1 > for all z G [0, oo), relation (2.13) gives a uniform estimate in (e,t) € 
(0, 1) X (0, T+(e)) of {dxfe,dxge) in L2{{0,T), L2 x L2) in dependence only of the initial condition 
(/o, (7o)- Indeed, on the one hand, since Inz < z — 1 for all z £ [0,oo), we have 

\2 



SO that 



zln z - z + 1 < z{z - 1) - {z - I) = {z - I 

R 



for all z > 0, 



£l{foe,90e) < j {{foe 



R,. 



(gos - ly ] dx<Ci 



(2.14) 



(2.15) 



for all e G (0, 1) by (2.4). On the other hand, owing to the Poincare-Wirtinger inequality and (2.9), 
we have 



ll/.(t)l|2 < 



fe{t)-^\\fe 



mm 



- < C iWdxfe 



+ 1) 



2 VL 

A similar bound being available for g,., we infer from (2.13), (2.15), and the nonnegativity of £1 
that, for T G (0,r+(e)), 



T 



{\\fe{tmi + \\9e{tWHi) dt < C / [l+Wdxfeimi + Wdxgemi) dt 

Jo 

< C [T + £i{f^e,goe)]<C2{T). (2.16) 

We next use this estimate to prove that the solution {fe,9e) of (2-2) is bounded in Loo{{0,T), x 
H^) for all T < T^{e). While the estimates were independent of e up to now, the next ones have a 
strong dependence upon e which explains the need of a regularisation of the original system. 
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Lemma 2.5. Given (e,T) £ (0,1) x (0,T+(e)), there exists a constant C{e,T) > such that the 
solution (fsjge) of (2.2) fulfills 

\\femm + \\ge{t)\\m<C{e,T) for all t e [0,T]. (2.17) 

Proof. We prove first the bound for f^. Given z G M, let q{z) := z'^/2. With this notation, the first 
equation of (2.2a) reads 

dtfs - (1 + R)dlq{fe) = Rd.iife - e)d,Ge). 
Multiplying this relation by dtq{fe) and integrating over (0,L), we get 



dtfedtqife) dx-{l + R) [ dlq{f,)dtq{fe) dx = R [ d,{{f, - e)d,G,)dtq{fe) dx. 

Jo Jo 

Using an integration by parts and Young's inequality, we come to the following inequality 
\\^/Te^tfe\\^ + ^Jt\\9-l(fe)g < \\\^/Te^tfe\\^ + feidAife " e)d,Gs)]' dx. 



Whence, we have shown that 



\\y/Tsdtfe\\l + (1 + R)±\\dM)\\l < R^ £ [fhdlG,)^ + feidM\d,Gsf] dx, (2.18) 
and the second term on the right-hand side of (2.18) may be estimated, in view of (2.8), by 

r fe{d,fef{d,G,fdx<- r f^{dJ,f{d.Gefdx = - f\dM))\d^G,f dx 
Jo ^ Jo £ Jo 



<\\\d,Ge\\Ud.q{fe)\\l 



Now, since G^ is the solution of G^ — s^d'^Ge = with homogeneous Neumann boundary conditions 
at X = 0, L, and g^ > e, the elliptic maximum principle guarantees that G^ > e. Hence, —e^d'^Gs < 
and therefore, for x G (0, L), 

/•X i-L 

-e'^dxGeix) < / gedx<\\gs\\i and e'^d^Gsix) < / gsdx<\\gs\\i, 
Jo Jx 

so that e^||9xGe||oo < lls'elli- In view of (2.9), we arrive at 

("^ fe{d,fe)\d,G,fdx<^\\g,\\l\\dM)\\l<C{eW^^^ (2.19) 
Jo ^ 

Concerning the first term on the right-hand side of (2.18), it follows from (2.8) and (2.12) that 
r flidlGefdx < - f ft{dlG,fdx < -\\q{fe)\\lo\\dlGe\\l < ^Mfe)\\lo\\dxge\\l 

Jo ^ Jo £ £ 

In order to estimate ||(?(/e)||ooi we choose S (0, L) such that ^/^(xe) = H/elli- By the fundamental 
theorem of calculus, we get 

q{fs){x) = q{fe){xe)+ jyxq{fe)dx < -i^||/,||2 + VZ||a,g(/,)||2 for ah x G (0,L). 
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Summarising, we obtain in view of (2.18) and (2.19), that 

J^\\d.q{fe)\\l < C{e)il + \\d,gs\\l){l + \\dMfe)\\l), (2.20) 
and from (2.16) and (2.20) 

||9.g(/,(t))||2<(l + ||a,g(/oe)||i)exp(^C(e)^*(l + ||a,5ei) ds^ < C{e,T), t e [0,T]. 
Since > e, we then have 

\\d^feml<C{e,T) foralHG [0,r]. 

Using (2.9) and the Poincare-Wirtinger inequahty, we finahy obtain that 

\\fe{t)\\m <C{e,T) for ant G [0,T]. 

Moreover, due to the symmetry of (2.2a), satisfies the same estimate as f^, and this completes 
our argument. □ 

Proof of Theorem 2.1. The proof is a direct consequence of Lemma 2.2, the lower bounds (2.8), and 
Lemma 2.5. □ 

We end this section by showing that, though £2 is not dissipated along the trajectories of the 
regularised system (2.2), a functional closely related to £2 is almost dissipated, with non-dissipative 
terms of order e. 

Lemma 2.6. For e G (0, 1) and T > 0, we have 

rT 

£2AfeiT),9e{T)) + / fs\{l + R)d.fe + Rd^.Ge\'' + RR^,ge\dAFe+ge)f dt 

Jo L ^ 

< S2Af0e,90e)+eC2 [ Qeit) dt, (2.21) 

Jo 

with Qe := WdxfeWl + ll^xSelll "^^^ 

2£2Me,9e):={l+R)\\fef2 + R\\gsf2 + R [^F.g, + Gefe) dx. (2.22) 

Jo 

Proof. We multiply the first equation of (2.2) by {1 + R)fs + RGg and integrate over (0, L) to obtain 

dtfe ((1 + R)fe + RGe) dx = - ^ f, ((1 + R)d^f, + Rd^G^f dx + Ii,, (2.23) 
Jo Jo 

with ^ 

/i,, := eR [ d^Ge ((1 + R)d^fe + Rd,G,) dx. 
Jo 

Thanks to Holder's inequality and (2.11), we have 

< eR{l + R) {Wd^G.hWd^feh + Wd^GeWl) < eG {\\dje\\l + 11^x5.111) • (2.24) 
Similarly, multiplying the second equation of (2.2) by R{Fi; + g^) and integrating over (0, L) give 

R / dtge {Fe + ge) dx = -RR^ / ge {d^Fe + d^gef dx + /2,e (2.25) 
Jo Jo 



A DEGENERATE PARABOLIC SYSTEM MODELLING FLOWS IN POROUS MEDIA 



11 



with 



l2,e ■■= eRR^, [ d^FA {Fe + ge) dx. 
Jo 

Using again Holder's inequality and (2.11), we obtain 

\h,e\ < eRR^ {\\d^F,\\l + Wd^QehWd^Feh) < eC {\\dM\l + \\dxge\\l) ■ (2.26) 



Observing that 

rL f-L 



[ dtfe ((1 + R)fe + RG,) dx + R [ dtge {F^ + ge) dx 

Jo Jo 

1 + ^ ^ - "2 ^ {dtFe - E^dldtFe) dx 

Jo 



2 dt''''' 



10 

-L 



+ ^Jthe\\l + Rj^ FeidtGe-e'dldtGe) dx 
= ((1 + R)\\f eWl + R\\9e\\l + 2R £ FeGe dx^ 

+ e^R [ {d,Gedtd,F, + d^F^dtd^Ge) dx 
Jo 

= ~(^il + R)\\fe\\l + R\\9ef2 + 2Rj^ {F,Ge + e^d,FAGe) dxy 
we sum (2.23) and (2.25), use (2.24) and (2.26) to obtain 

((1 + R)\\fe\\l + RheWl + 2R £ {F,Ge + E^d^FAG,) dx^ 

< - r fe ((1 + R)dJs + Rd^Gef dx - RR^ ge (d.Fe + d.gef dx + eC2 Qe- 
Jo Jo 



Since 



dx 



2 [ [FeGe + e'd^Fed^Ge) dx = [ (2F,G, - e^GedlFe - e^FedlGe) 
Jo Jo 

= {Gefe + Fege) dx, 

Jo 

the claimed inequality follows from the above two identities after integration with respect to time. 

□ 

3. Weak solutions 

Given T G (0, oo], we let Qt := (0,T) x (0,L). Furthermore, given e G (0,1), we let {fe,ge) be 
the global strong solution of the regularised problem (2.2) constructed in Theorem 2.1. We shall 
prove that {fe,ge) converges, in appropriate function spaces over Qt, towards a pair of functions 
{f,g) which turns out to be a weak solution of (1-4) in the sense of Theorem 1.1. 
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Recall that, by (2.8), (2.4), (2.9), and (2.13), {fe-,9e) satisfies the following estimates 

(o) /e > e, 9e>^ on Qoo, 

{h) Wheh + Iboelb < II/0II2 + II90II2 + 2VL, 

ae{t)\\i = llffolll + ^L, 
R 



(c) ||/e(t)||i = ||/o||i+eL, 

{d) £i{fe{t),ge{t)) + 
for t > 0. Using (3.1), we show that: 



(3.1) 



s\\2 



+ 



1 + 2i?' 



Lemma 3.1 (Uniform estimates). Let h £ {f,g,F,G}. There exists a positive constant C4(T) such 
that, for all (e,T) G (0, 1) x (0, 00), we have 

+ \\heml) dt<C^{T), (3.2) 



(0 







2 



{^^) 



T 



\dtK 



ll^^ydt<C,{T). 



(3.3) 



Proof. The estimate for in L2{0,T, H^) is obtained from the energy estimate (2.13), by taking 
also into account relations (2.9), (2.10), (2.11), (2.15), the nonnegativity of £1, and the Poincare- 
Wirtinger inequality as in the proof of (2.16). In order to prove the second estimate of (3.2), we 
note that, since is continuously embedded in Lqo and (/g) is bounded in L2{0,T, H^), (/g) is 
uniformly bounded with respect to e in Loo(0,oo;Li) n L2(0,T;Loo) for all T > 0. The claimed 
L3— bound then follows from the inequality WfeW^ < II /e II ^ II /e 111- Next, an obvious consequence of 
the definition of is that ||-Fe||p < ||/e||p foi" P ^ [l;Oo], from which we deduce the expected bound 
in L3 for (-Fe). A similar argument shows that (g^) and (Gg) satisfy the second estimate in (3.2). 

In order to prove (ii), consider first h € {f,g}. From (3.2) and Holder's inequality we obtain that 
ifsdxfe), {{f£-£)dxGe), {{ge-e)dxFs), and {ged^ge) are uniformly bounded in Lg/sl^r)- Therefore, 
the equations of (2.2a) may be written in the form S^/ig = dxH^, for some function which is 
uniformly bounded in Lq/^{Qt) and satisfies homogeneous Dirichlet conditions H^{0) = Hl'{L) = 0. 
Given 



€ Wq, we have 



\{dthe\(t>)i 



Consequently, 



\dthe 



(pdxH^ dx 



H^dx(l)dx 



< 



6/5 



for t £ (o,r). 



and the families (/e), (gs) are both uniformly bounded in Lq/^{0,T; {WqY). 

Finally, we have ((1 — e^d^)p\q)L2 — (p\{^ ~ ^'^(^x)q)l2 for all p,q £ H^, and choosing p :- 

(1 - £2, - ' - " ' 



'^dtfe, q = {1 — e^d^) ^(p with (f) G Wq, we have 

{dtFe\cP)L, =((1 - e^dlr'dtfe\<P)L, = {dtfe\{l 



--{dxHf\{l-s'dl 



)L2 



-{Hf\il-e'dl)-'dx 



)L2- 



A DEGENERATE PARABOLIC SYSTEM MODELLING FLOWS IN POROUS MEDIA 13 
Since (1 — e^9^)~^ is a contraction in Lq, we obtain that 

mFe\^)L,\ < \\Hf\\,/,\\{l-e'dlr'd,ne < Il^/ll6/5ll</'llw^i, 

and the assertion (ii), when h = F, follows at once. Invoking a similar argument for (G^), we 
complete the proof. □ 

This lemma enables us to use a result from [9] and show that (/e), (fi'e), (^e), and (Gg) are 
relatively compact in L2(0, T; G°([0, L])), provided that a G (0,1/2). This will allow us to identify 
a limit point for each of these sequences, and find in this way a candidate for solving (1.4). Indeed, 
we have: 

Lemma 3.2. Given h £ {f,g,F,G}, and a S (0,1/2), there exists a subsequence [hf,^) o/ (/i^) 
which converges strongly in L2(0, T; G"([0, L])). 

Proof. Invoking the Rellich-Kondrachov theorem, we have the following sequence of embeddings 

H' ^C''{[Q,L])^{Wl)\ a<l/2, 

with compact embedding G°([0, L]). Furthermore, in view of Lemma 3.1 (i), the family {he) 

is uniformly bounded in L2{0,T; H^), while, by Lemma 3.1 (ii), {dthe) is uniformly bounded in 
Li(0, T; (VFg^)'). Whence, the assumptions of [9, Corollary 4] are all fulfilled, and we conclude that 
(he) is relatively compact in L2(0, T; G°([0, L])). □ 



3.1. Construction of weak solutions. Using the uniform estimates deduced at the beginning of 
this section, we now establish the existence of a weak solution of (1.4). Owing to Lemma 3.2, 
there are f,g,F,G £ L2(0, T; G"([0, L])) such that, for a £ (0,1/2), 

fe.^f, ge.^g, Fe^^F, Ge.^G inL2(0,r;C"([0,L])). (3.4) 

Furthermore, by Lemma 3.1 (i), the subsequences {dxfe^), {dxge^), (dxFe^), and [d^Gs^) are uni- 
formly bounded in the Hilbert space L2{Qt)- Hence, we may extract further subsequences (denoted 
again by (/^J, {g^^), (F^J, and (G^J) which converge weakly: 

dxfek 9xf, dxgek 9xg, d^F^^ d^F, d^Ge^ d^G in L2{Qt)- (3.5) 

In fact, we have that 

f = F and g = G a.e. in Qt- (3-6) 

Indeed, (3.6) follows by multiplying the relation F^^. — el^^F^j. = /^^ by a test function in H^, 
integrating by parts, and letting then /c — ?• oo with the help of (3.4) and (3.5). In view of (3.4)- 
(3.6), we then have 

fsk^xfek fdxf, fe.dxGs^ ^ fd^g in Li{Qt), 

(3.7) 

gskdxFe^ gdxf, ge^dxgek gdxQ m Li{Qt). 
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Using the fact that {fe,ge) are strong solutions of (2.2), we obtain by integration with respect to 
space and time that 

f-L f-L 



[ fe.mdx- [ foe,lPdx=- [ [{l + R)f,^d^f,^+R{f,^-ek)d:,Ge,] d^ijdxdt, 

Jo Jo JQt 

/ ge^iT)ipdx- goekipdx=-R^ [{ge^ - ek)dxFs^ + ge^d^^gek] ^xi^ dx dt, 

Jo Jo JQt 



for all T > 0, e G (0, 1), and S TV^. Since 

/oefc fo and go^^ go in L2 (3.9) 

by classical arguments, we may pass to the limit as A; —t- 00 in (3.8) and use (3.4), (3.5), (3.7), 
and (3.9) to conclude that (/,<?) is a weak solution of (2.2) in the sense of Theorem 1.1. The fact 
that (/, g) can be defined globally follows by using a standard Cantor's diagonal argument (using a 
sequence r„ 00). 

3.2. Energy estimates for weak solutions. Letting e — )• in the relation (3.1) (c), we find in 
view of (3.4), that 

ll/(i)lli = ll/olli, Il5(t)l|i = Ibolli, tG(0,oo). 
We show now that the weak solution found above satisfies the energy estimate 



£i{f{T),g{T)) + Q|5^/|2 + _JL_|a^5|2^ dxdt<£r{h,go) 
for T E (0,00). Recall that, by Lemma 2.4, we have 



(3.10) 



fl(/,,(T),5,,(r))+ j {l\d^f^f + -^—\d,g^f^ dxdt<£,{foe„gOe,) (3.11) 



for all A: G N. On the one hand, note that (3.4) and Fatou's lemma ensure that 
while (3.5) implies 



£i{f{T),g{T)) < limmi£i{f,^{T),g,^{T)) for T G (0,oo), (3.12) 



/ \dx f l"^ dx dt <limmi / {dxfej'^dxdt, 
JQt JQt 

j l^i^^p (ix < liminf / \dxgej\^ dx dt. 
Jqt Jqt 



(3.13) 



We still have to pass to the limit in the right-hand side of (3.11). By (3.9), we may assume that 
(/osfc) and (goek) converge almost everywhere towards /o and go, respectively. Furthermore, since 
< xlnx — X + 1 < 1 + 2x3/2 fo^ 

X > 0, we have 

/ \fos,lnfoe,-foe, + l\dx<\E\+2 [ f^^^^dx<\E\+2\E\'/^\\fos,\\T <C\E\'^^ 

J E J E 

for all A; G N and all measurable subsets E of (0, L), meaning that the family (/oej. In /oe^, — /oej, + 1) 
is uniformly integrable. Clearly, the same is true also for {gock doe^ ~ 9oek + infer then from 
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Vitali's convergence theorem, cf. [5, Theorem 2.24], that the hmit of the right-hand side of (3.11) 
exists and 

lim £i{foe^^,goef,) =£i{fo,9o)- 

Whence, passing to the Umit in (3.11), we obtain in view of (3.12) and (3.13) the desired estimate 
(3.10). 

Finally, we show that weak solutions of (1.4) satisfy 



£2{f{T),g{T))+ f{{l + R)dxf + RdxgY + RR^g{dJ + dxgY dxdt < £2{fo,go) (3.14) 
■JQt ^ ^ 

for T G (0, oo). In virtue of (3.4), (3.5), and (3.6) we have 

\fh'k^\fl^ dxfek^dxf, and dxGe,, ^ d^g in L2{Qt) 
which imphes that ^fU'uQxfsk ^ VJdxf and ^/J^dxCs,^ VJdxg in Li{Qt). Consequently, 

^,{{1 + R)dxfe, + RdxGe,) - 77((1 + R)dxf + Rdxg) in Li{Qt). 
and, by a similar argument, 

y/g^{dxFe^ + dxgsk) Vaidxf + dxg) in Li(Qr)- 

Now, owing to (3.2), the sequence {Qsk) defined in Lemma 2.6 is bounded in L2((0, T)) and we 
then infer from Lemma 2.6 that both (^^((1 + R)dxf£^ + RdxGs^)) and {^/g^{dxFs^ + dxge^)) 
are bounded in L2{Qt)- The previous weak convergences in Li{Qt) may then be improved to weak 
convergence in L2{Qt) (upon extracting a further subsequence if necessary) and we can then pass 
to the limit in (2.21) to conclude that (3.14) holds true, using weak lower semicontinuity arguments 
in the left-hand side and the property EkQe^. — )• in L2((0, T)) in the right-hand side. 

3.3. Exponential convergence towards equilibria. In this last part of the paper we prove our 
second main result. Theorem 1.3. The proof is based on the interplay between estimates for the two 
energy functionals £i and £2, with the specification that we use £1 to estimate the time derivative of 
the stronger energy functional £2, and obtain exponential decay of weak solutions in the L2— norm. 
Recall from (3.1) (c) that 



Ak :-- 



Olll 



L L 
for all G N and t E [0, 00). Introducing 



+ e, = M^ and B,:=Ml + ,, 



/..In 



R 
Rii 



9e, In 



Bk 



9e, + Bk 



dx 



+ l£ {ife, - Akf + R [ifs, - Akf + {ge, - Bkf + (ge, - Bk){Fe, - Au) 
+{f,^-Au){Ge,-Bu)]}dx, 
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we infer from (3.1) (c) and the proofs of Lemma 2.4 and 2.6 that 



dt 



R 

dtfsk In/efe + -^dtQe^ In 



dx 



^\jt[ ^^^^ ^^^'^ + ^ [^'^ + ^^^^^^ + ^ 

5eJli + efcC'3 (ll^x/ejli + \\dx9ek\\l) 



1 

2' 
1 

3' 



<-^l|5x./.jii 



/2 



1 + 2i?' 
R . 



%9ek Hi 



2 + 2i?' 

provided k is large enough. Using the Poincare-Wirtinger inequahty, we find a positive constant 
C5 such that 

d-^k ^ ^ ^iij, ^i|2,|| Dl|2^ 



<-C5 (ll/e, -Afc||i + ||g,, -Sfclli) 



(3.15) 



for large k and ah t € (0, 00). We show now that the right-hand side of (3.15) can be bounded by 
—ojJ-k for some small positive number uj. Indeed, arguing as in Lemma 2.3, we find that 

[{ge, - B^m^ - Ak) + (/e, - Ak){G,, - Bk)] dx 

Jo 



< 



\9ek 



BkhWFe, - Akh + Wfe, - AkhWG. 



B^ 



< 2\\g,^ - BkhWfe, - Akh < Wfe, - AkWl + \\9e. 
Recalling (2.14), we end up with 

r-L / 



fc 2 



BkWl- 



(3.16) 



\\f,^-Ak\\i + \\ge,-B. 



k\\2 



Al 







Ak 



> 



>min <j Ak, 

k 



Ak 

RfiBk 
R 



Ak 

L 



Ak 



+ Bk 



At[^ln[^\-^ + l\+BU^ln[^\-^ + l 



9ek _ 

Bh 



9eu 



Bk \Bk 



dx 



9ek 



Bk 



dx 







A. 



R 



/..In - fe,+A, + _ ^,^^ln -9e,+B, 



Rij 



Bu 



Combining (3.15), (3.16), and (3.17), we conclude that if ||/o||i > and ||(7o||i > 0, then 

dFk 



dx. 
(3.17) 



dt 



it) < -ujFk{t) 



for some positive constant uj and k sufficiently large. Whence, 

\\h^{t) - AkU, + \\9e,{t) - BkU, < Ce-^\ 

which yields, for k — )• 00, the desired estimate by (3.4) and (3.6), as stated in Theorem 1.3. 

If /o = [resp. gQ = 0], then / = [resp. g = 0], while g [resp. /] is a weak solution of the 
one-dimensional porous medium equation and converges therefore even in the Loo— norm to flat equi- 
libria (if /o = 0, then uniqueness of solutions to (2.6) implies that /e = e and dtge = Rii{g£dx9e))^ 
cf. [10, Theorem 20.16]. Convergence in this stronger norm is due to the fact that comparison 
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methods may be used for the one-dimensional porous media equation, while for our system they 
fail because of the structure of the system. 
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